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Abstract 

SeCQC is an open-source program code which implements a Nu- 
merical Search for the classical Capacity of Quantum Channels 
(SeCQC) by using an iterative method. Given a quantum chan- 
nel, SeCQC finds the statistical operators and POVM outcomes 
that maximize the accessible information, and thus determines 
the classical capacity of the quantum channel. The optimization 
procedure is realized by using a steepest-ascent method that fol- 
lows the gradient in the POVM space, and also uses conjugate 
gradients for speed-up. 



This manual is for version 1.0. 
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1 License Agreement 

SeCQC is an open-source program that, given a quantum channel, imple- 
ments a Numerical Search for the classical Capacity of Quantum Channels. 
It is a derivative of the open-source program code SOMIM (see Ref. [1]). 
Copyright © 2010 J.W. Shang, K.L. Lee and B.-G. Englert. 

SeCQC is a free software: You can redistribute it and/or modify it under 
the terms of the GNU General Public License Version 3 as published by the 
Free Software Foundation. 

SeCQC is distributed in the hope that it will be useful, but WITH- 
OUT ANY WARRANTY; without even the implied warranty of FITNESS 
or MERCHANTABILITY FOR PARTICULAR PURPOSE. See the GNU 
General Public License at http://www.gnu.org/licenses/ for details. 

2 What can SeCQC be used for? 

Kraus representation of a Quantum Channel: We use a set of Kraus 
operators K m (m = 1, N) to represent the quantum channel, such that 

JV 

£AtA m = l. (1) 

m=l 

Consider the following quantum communication scenario shown in Fig. 1. 
Alice sends a set of quantum states £ = {pj \ j = 1, 2, J} with pj > to 
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Bob through a quantum channel M, such that 

J 

p = Y.P3 withtr{p} = l. (2) 

3=1 

The states after passing through the channel are given by 

TV 

P ;. = £ K mP] Kl = Mp 3 . (3) 

m=l 

Bob performs a generalized measurement, specified by a positive-operator- 
valued measure (POVM), on the state he receives. The POVM with outcomes 
Vi k {k = 1, 2, K) decomposes the identity, 

K 

]Tn fc = i withn fe >o. (4) 

k=l 

Then, the joint probability to receive the j'th state and get the kth outcome 
is 

Pjk = tripjUk} , J2p jk = 1 • (5) 
Bob's figure of merit is the mutual information 

^;n) = EE^iog 2 ^r- (6) 

j=l k=l Pj-P-k 

where pj. and p.k are the marginal probabilities, 

Pi- = Y,Pik = tr{pj} , p. k = YsPik = E tT {p( K l u kK m )} . (7) 

k j m 
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As stated, the pjs are normalized such that their traces equal the probabilities 
of receiving them. 

Classical Capacity of Quantum Channels: Generally, the classical 
capacity of a quantum channel can be defined as the maximum accessible 
information with respect to both statistical operators and POVM outcomes, 

C = max max I(E ; II) . (8) 
s n 

Given a certain quantum channel, SeCQC finds the statistical operators as 
well as POVM outcomes that maximize the accessible information (AI), and 
thus determines the classical capacity of the quantum channel. 

The calculation is performed using a combination of the steepest-ascent 
method (see Ref. [2] and Section 11.5 in Ref. [3]) and the conjugate-gradients 
(CG) method [4]. The percentage chance to calculate with one method or the 
other can be specified by the user (see Section 4 below). The implementation 
in SeCQC also makes use of the golden-section search method. 

3 Download and Compile 

The complete set of files, including this manual, are available at the SeCQC 
site: http : //www . quantumlah . org/publi cat ions/ sof tware/SeCQC/. Down- 
load http : //www . quantumlah . org/publications/ sof tware/SeCQC/ all . tar . 
gz, if you want to have the complete collection of files. Just this manual is 
fetched from http : //www . quantumlah . org/publications/sof tware/SeCQC/ 
Manual.pdf. The Windows executable file for SeCQC can be downloaded 
from http : //www. quantumlah. org/publications/sof tware/SeCQC/ secqc . 
tar.gz. If you intend to modify the code, you can download the source files 
from http : //www. quantumlah. org/publications/sof tware/SeCQC/ source . 
tar . gz. 

The program is written in C++ and the graphic user interface (GUI) 
is implemented using wxWidgets (http://www.wxwidgets.org/). Here are 
the instructions for compiling SeCQC: 

1. Install wxWidgets from http://www.wxwidgets.org/downloads/. 

2. If you are working in Windows, you need to install MinGW (http:// 
www.mingw. org/download . shtml) and MYSY (http : //www.mingw. org/ 
msys . shtml) as well. 

3. When wxWidgets and MinGW are configured, you can compile SeCQC 
by executing g-| — |- MI.cpp 'wx-config —libs' 'wx-config — cxxflags' 
-o YourProgramName in MSYS shell. 
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4. If you face problems running the program in a Linux environment, try 
export LD .LIBRARY _PATH= /usr /local /lib. 

5. The executable file is compiled under Windows XP Service Pack 3, with 
wxWidgets 2.8.10, MinGW 5.1.6 and MSYS 1.0.11. 

4 How to use the program 

The GUI of SeCQC is shown in Fig. 2. In the first box labeled as "Parameter 
Settings", J is the number of statistical operators pj. The current maximum 
possible value is J = 30. Parameter K is the initial number of POVM 
outcomes, with the largest possible value being K = 30. The third and fourth 
fields are the number C of Quantum Channels to be inserted and the number 
M of Kraus operators for each channel respectively. The next two fields 
are the dimension D of Kraus operators and the dimension N of statistical 
operators respectively with 30 being their highest possible value. All the 
maximum values mentioned above can be changed by modifying the source 
code. The seventh field is the percentage chance to use the steepest-ascent 
method to perform maximization in an iteration; this parameter controls the 
relative frequency of using the direct or the conjugate gradient. The eighth 
field gives the tolerance in the accessible information, the stopping criteria 
for the computation; the calculation stops when the difference in accessible 
information between the current iteration and the previous iteration is less 
than half of the sum multiplied by the tolerance plus the machine_epsilon 
e m (also termed as the machine accuracy, typical value for double precision 
is around 1.6 x 10~ 16 ), i.e. when 2.0 x (current — previous) < tolerance x 
(current + previous) + e m . The ninth field is the name of the output file. 
By default, the output file will be located at hard disk C. You can change 
the output directory by clicking the ellipsis button and choose your 

preferred location. 

The next three boxes display the input Kraus operators {A m } m=1 . jM , the 
optimal statistical operators {Pj}j=i,...,j and the calculated optimal POVM 
outcomes {Hk}k=i,...,K- The spin buttons are used to switch between the 
various K m s/ pjs/UkS, while the small box beside the spin button is used to 
choose to display the real or imaginary part of the chosen K m s/ pjs/UkS. 

The maximum accessible information for the given channel will be dis- 
played in the last box after the "Calculate MI" button is pressed. All values 
will be reset to default when the "Reset" button is pressed. 

Important note: The matrices for the K m s must have the correct di- 
mension; they must satisfy the condition, i.e. J^mK^Km = 1. 
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jfc SeCQC Version 1,0 



File Cak About 
input, bet | 
Parameter Settings 



No. of Quantum States 3 -■ 
No, of Quantum Channels C : 
Dimension of p N = 
Percentage chance using DG = 
Output file = 



output, bet 



T| No . of POVM outcomes K = 
~~{ No. of Kraus Operators M = 
jx 2 Dimension of Cap D = 

~7"j Tolerance in MI - 

Save output file to 



Input Matrix Representation of Kraus Operators 
Matrix representation of Cap |"l 



Display Real or Imaginary — 
(* Real P Imaginary 



0.0000 


D.8660 


0,8660 


0.0000 



Display Optimal Quantum States 



Display optimal density operator p |T" 



Display Real or Imaginary — 
[? Real P Imaginary 



0,2500 


-0.2500 


-0,2500 


0.2500 



Display Optimal POVM outcomes 



Display the optimal n [T~ 



Display Real or Imaginary — 
(* Real C Imaginary 



0.5000 
0,5000 



0.5000 
0.5000 



Display the Maximal Mutual Information 



Maximal Mutual Information I = 1 1.000000 
Calculate MI Reset 



Figure 2: The graphical user interface (GUI) of the program. 
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5 How to import data 

Data can be imported into SeCQC using a text file that is possibly generated 
by another program. An example is shown in Fig. 3. When importing, the 
numbers after the equal signs will be read into the program. The first line is 
the dimension D of the Kraus operators. The second line gives the dimension 
N of the statistical operators. The third line gives the number J of the 
statistical operators and the fourth line gives the number of outcomes K of 
POVMs that the program should start calculating with. And the last two 
lines give the number M of the Kraus operators for each channel and the 
number C of Quantum Channels respectively. 



W input.txt - Notepad 




_|n|x| 


File Edit Format View Help 


D= 2 








N= 2 








3= 2 








K= 5 








M= 2 








C= 1 








Capl = {{0.5, 0}, {0, 


5}} 






Cap2 = {{0, 0.8660254}, 


{0.8660254, 0}} 


















J 2i 





Figure 3: Example of an import file. 



The subsequent lines give the input matrices for the Kraus operators. 
Each line will give only one operator. For an operator represented in matrix 
form as 

0.1 0.3 + 0.5i\ (] 
0.3-0.5z 0.6 J' W 

the input data should be formatted as {{0.1, 0.3 + 0.51}, {0.3 — 0.51, 0.6}}. 

Complex numbers are entered as RealPart + ImaginaryPart I, as illus- 
trated by —3.1 — 4.51. Please note that the complex unit % must be entered 
in upper case I and it must be at the end of the entry. 

6 Meaning of output data 

A typical output file looks like Fig. 4. The first eight lines give the follow- 
ing information: the number J of statistical operators, the initial number 
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File Edit Format View Help 

Number of possible Quantum States J = 2 
Number of POVM outcomes K = 5 
Dimension of statistical operators N = 2 
Number of Kraus Operators Cap_m = 2 
Number of Channels Num_channels = 1 
Dimension of K_m = 2 

Percentage chance to calculate with direct gradient = 0/100 
Tolerance in Mutual Information = le-006 
Starting seed for number generator = 1285115359 



counte r 




1, 


new_HI = 


0.169736 


counter 




2, 


new_MI = 


0.188696 


counte r 




3, 


new_MI = 


0.18872 


counte r 




4, 


new_MI = 


0.188725 


counte r 




5, 


new_MI = 


0.188733 


counte r 




6, 


new__MI = 


0.188847 


counte r 




7, 


new__MI = 


0.189853 


counte r 




8, 


new__MI = 


0.197004 


counte r 




9, 


new__MI = 


0.243591 


counte r 




10, 


new__HI 


= 0.438612 


counter 




11, 


new__HI 


= 0.71935 


counter 




12, 


new__HI 


= 0.889554 


counter 




13, 


new__HI 


= 0.963343 


counter 




14, 


new__HI 


= 0.98875 


counter 




15, 


new__HI 


= 0.996663 


counter 




16, 


new__MI 


= 0.999022 


counter 




17, 


new__MI 


= 0.999709 


counter 




18, 


new__MI 


= 0.999907 


counter 




19, 


new__MI 


= 0.999968 


counter 




20, 


new__MI 


= 0.999987 


counter 




21, 


new__MI 


= 0.999996 


counter 




22, 


new__MI 


= 0.999999 


counter 




23, 


new__MI 


= 1 



The Kraus Operators of the Capacitor are: 
Capl={{0. 5,0}, {0,0. 5}} 
Cap2={{0, 0. 866025} , {0. 86602 5 , 0}} 

The initial input states are: 
pl={{0. 2 50237,0}, {0,0.28 5306}} 
p2={{0. 246993,0}, {0,0.217464}} 

After optimization, the input states are: 

pl={{0. 249934 -1. 58681e -0181, 0.249995 -3. 71209e -0061}, {0.249995+3. 71209e -0061, 0.2 
p2={{0. 2 50065+1. 5874 5e -0181, -0.2 50005+4. 06794e -0061}, {-0.250005 -4. 06794e -0061,0 

The initial POVM outcomes generated are: 

nl={{0. 49995 5,0. 499901+1. 29101e -0051}, {0.499901-1. 29101e -0051, 0.499846}} 
n2={{0. 499944, -0. 5-1. 24069e -0051}, {-0. 5+1. 24069e-005I, 0. 500056}} 
n3={{2.7218e-006,2.62809e-006-2.86621e-008l}, {2 . 62809e-006+2 . 86621e-008I, 2 . 5379 
n4={{4.95261e-005,4.95078e-005-1.16994e-007l}, {4 . 95078e -005+1. 16994e-007I ,4 . 948 
n5={{4.82 3 58e-005,4.73167e-005-3. 575 32e-007l} , {4 . 73167e -005+3 . 575 32e -0071, 4. 641 

After eliminating the extra POVM outcomes: 

nl={{0. 500056, 0.5+1. 24069e -0051}, {0. 5-1. 24069e-005l, 0.499944}} 
n2 = {{0. 499944, -0. 5-1. 24069e -0051}, {-0. 5+1. 24069e-005l, 0. 500056}} 




Figure 4: Example of output file. 
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K of POVM outcomes, the dimension N of the statistical operators, the 
number M of Kraus operators for each channel, the number C of Quantum 
Channels inserted, the dimension D of the Kraus operators, the tolerance 
in the calculated mutual information, and the seed for the random number 
generator. 

The next block of lines gives the mutual information at the end of each 
iteration. In the example shown in Fig. 4, altogether 26 iterations have been 
performed with the final accessible information being AI = 1.000 exactly. 

The subsequent five blocks of lines give the M Kraus operators for each 
channel, the initial as well as the optimized J statistical operators and the K 
outcomes of the optimal POVMs that correspond to the accessible informa- 
tion calculated in the final round of iteration. Each Kraus oeprator/statistical 
operator/POVM outcome is given in a single line in matrix form, as explained 
in Section 5. 

Among the K outcomes, if any two outcomes, say IT^ and ITfc 2 , give 
equivalent probabilities, i.e. Pjk 1 P-k 2 — P-k 1 Pjk 2 f° r & U j, then these two 
POVM outcomes are replaced with one new POVM outcome, + n fe2 , 
such that the new optimal POVM contains only K — 1 outcomes. The last 
block of data in the output file gives the POVM after this elimination process, 
i.e. the POVM is the optimal POVM with the least number of outcomes. 

Caution: As is the case for all steepest-ascent methods, there is the 
possibility of convergence towards a local, rather than a global, maximum. 
There is no absolute protection against this danger, but in practice one can 
fight it efficiently by running the program many times for comparison, with 
different seeds. It also helps to start with a rather large K value. 

7 Contact information 

Please send your comments, suggestions, or bug reports to the following email 
account: secqc@quantumlah . org 
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